Novel universality class of absorbing transitions with continuously varying critical 

exponents 



O 

o 

o 

Q 

o 



o3 



13 
G 
O 

o 



> 
in 

On 
O 



i 

G 

O 

o 



X 



Jae Dong Noh 

Theoretische Physik, Universitdt des Saarlandes, 66041 Saarbriicken, Germany 

Hyunggyu Park 

School of Physics, Korea Institute for Advanced Study, Seoul 130-012, Korea 
(Dated: February 1, 2008) 

The well-established universality classes of absorbing critical phenomena are directed percolation 
(DP) and directed Ising (DI) classes. Recently, the pair contact process with diffusion (PCPD) 
has been investigated extensively and claimed to exhibit a new type of critical phenomena distinct 
from both DP and DI classes. Noticing that the PCPD possesses a long-term memory effect, we 
introduce a generalized version of the PCPD (GPCPD) with a parameter controlling the memory 
effect. The GPCPD connects the DP fixed point to the PCPD point continuously. Monte Carlo 
simulations show that the GPCPD displays novel type critical phenomena which are characterized 
by continuously varying critical exponents. The same critical behaviors are also observed in models 
where two species of particles are coupled cyclically. We suggest that the long-term memory may 
serve as a marginal perturbation to the ordinary DP fixed point. 

PACS numbers: 64.60.-i, 05.70.Ln, 64.60.Ht 



I. INTRODUCTION 

A noncquilibrium system with trapped (absorbing) 
states may display a so-called absorbing phase transi- 
tion between an inactive and an active phase 0, |^. A 
state which has a zero transition probability into any 
other state is called the absorbing state. A system in the 
inactive phase always evolves into the absorbing state 
and stays there forever. On the other hand, a system 
in the active phase may not be trapped in the absorbing 
state with a finite probability. There has been growing 
interest in the critical behaviors of the absorbing phase 
transitions since a wide range of phenomena such as epi- 
demic spreading, catalytic chemical reactions, and sur- 
face roughening, display absorbing transitions pi || . 

Besides their wide applications, absorbing critical phe- 
nomena have been the focus of a number of theoreti- 
cal works, since they are categorized into a few uni- 
versality classes characterized by the symmetry between 
the absorbing states and/or the conservation in dynam- 
ics H H Hi [§■ Criticality of each universality class is 
described by three independent critical exponents; f3 for 
the order parameter, v± for the correlation length, and 
u» for the relaxation time. For systems which are free 
from quenched disorder and evolve only through short- 
range processes, the directed percolation (DP) and the 
directed Ising (DI) (or parity conserving) classes are well- 
established ones. 

The DP class involves typically a single absorbing state 
without any kind of conservation in dynamics ^| . The 
contact process, a model for epidemic spreading, is a pro- 
totypical example of the DP class ||. In this model, 
individuals on a lattice arc either infected or healthy. In- 
fected ones may be healed spontaneously or infect healthy 
neighbors. There is a single absorbing state where all in- 
dividuals are healthy. Varying the relative rates between 



infection and healing processes, one can find a phase tran- 
sition from the absorbing phase into the active phase. 

Its stationary and dynamic critical behaviors are char- 
acterized by spatiotemporal cluster patterns of infected 
individuals. These clusters can be mapped to the di- 
rected percolation clusters 0, when the temporal di- 
rection is set to be the preferred direction of DP clus- 
ters. So all critical exponents (fractal dimensions) take 
the same values as the corresponding DP critical expo- 
nents. Most of systems with absorbing states belong 
to this DP class, e.g. the Domany-Kinzel cellular au- 
tomaton poL the Ziff-Gulari-Barshad model for a sur- 
face catalytic reaction the branching-annihilating 
random walks with an odd number of offspring jl^ ] , and 
the pair contact process (PCP) Unlike others, the 
PCP has infinitely many absorbing states, which leads to 
controversial transient behaviors, i.e., nonuniversal scal- 
ing H H, [H II versus absence of scaling 0, ||. But 
its stationary critical behavior still belongs to the DP 
class at low dimensions fl9f| . 

The DI class includes systems with two equivalent ab- 
sorbing states with Ising-like Z2 symmetry or equiva- 
lently in one dimension (ID) a single absorbing state 
with parity conservation in the domain wall language. 
The noncquilibrium kinetic Ising (NKI) model with com- 
bined zero-temperature spin-flip dynamics and infinite- 
temperature spin-exchange dynamics is an example of 
the DI systems f2C|| . In this model, only Ising spins near 
domain walls can flip, so that the two states with all spins 
up or down are absorbing. These two absorbing states 
are probabilistically equivalent. In terms of the domain 
wall, there is a single absorbing state (vacuum) with par- 
ity conservation in the number of domain walls, since spin 
flips change it only in pairs. Other examples in the DI 
class includes the interacting monomer-dimer model , 
the branching-annihilating random walks with an even 
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number of offspring 
the contact process 



and generalized versions of 



There also exist models in the 
DI universality class that have infinitely many absorbing 
states HJ, . They also display non- universal scal- 

ing behaviors in the transient regime p5| p(| . 

There is an infinite dynamic barrier between two ab- 
sorbing states of DI systems, which is similar to the free 
energy barrier between two ground states in the ordered 
phase of the equilibrium Ising system [Q . A state near 
one absorbing state cannot evolve into a state near the 
other absorbing state by a finite number of successive lo- 
cal changes. In other words, a frustration (domain wall) 
in a configuration generated by pasting two absorbing 
configurations cannot disappear within a finite number 
of time steps. For example, a spin state in the ID NKI 
model with all spins up in one semi-infinite lattice and all 
spins down in the other semi-infinite lattice, never relaxes 
to the absorbing state. 

The concept of the infinite dynamic barrier is very 
useful to understand the critical behavior of systems 
with infini tely many absorbing states. For instance, 
the PCP |p"3f and the modified interacting monomer- 
dimer (IMD-IMA) model p4| have infinitely many ab- 
sorbing states. In the PCP, a frustration between any 
of the absorbing states can disappear locally, so the ab- 
sorbing transition falls into the DP class. However, in 
the IMD-IMA, the infinite dynamic barrier separates the 
absorbing states into two equivalent groups of absorbing 
states, which results in the Dl-type critical behavior 

A few exceptional cases have been reported: Levy-type 
long-range flights are relevant to the absorbing critical 
phenomena, which lead to continuous variation of critical 
exponents (27j]. Multi-species particle reaction-diffusion 
systems also show non-DP and non-DI critical behaviors, 
where an inter-species hardcore interaction plays an im- 
portant role pq] . A lattice gas model with a global con- 
servation in the particle number shows a novel type crit- 
icality related to self-organized critical systems |^] . The 
nonequilibrium g-state Potts models in higher dimensions 
show interesting critical behaviors |30[ . A quenched ran- 
domness also leads to novel absorbing critical phenom- 
ena 

Recently, Howard and Tauber introduced a modified 
PCP (called as PCPD) model, which allows single par- 
ticle diffusion f|, ||. They studied the PCPD in the 
context of bosonic field theory and showed that the field 
theory is non-renormalizable and the absorbing transi- 
tion does not belong to the DP class. A fermionic version 
was first studied by Carlon, Henkel, and Schollwock 
in ID, which raises continuing debate on the universality 
class [35, 36, 3^, [3j| |3j| Hereafter, we only focus 
on the ID systems. 

The PCPD is defined on a lattice, each site of which is 
either occupied by a particle (X) or empty (0). Dynamic 
rules are given as follows. A nearest neighbor pair of par- 
ticles can either annihilate (XX — > 00) with probability 
p(l — d) or branch one offspring (XX — > XXX) to one 
of the neighboring sites with probability (1 — p)(l — d). 



A particle can hop to its neighboring site with probabil- 
ity d. Branching and hopping attempts are rejected if a 
particle would land on the top of another particle. 

When the diffusion is not allowed (d = 0), it reduces to 
the ordinary PCP with infinitely many absorbing states. 
Any state with only isolated particles is an absorbing 
state. But there is no infinite dynamic barrier between 
them and the transition belongs to the DP class pq |. 

At nonzero 0?, the PCPD has only two absorbing states; 
a vacuum state and a state with a single diffusing particle. 
In fact, the latter forms an absorbing subspace consist- 
ing of L configurations with system size L in ID. Once 
the system evolves into this absorbing subspace, it drifts 
freely inside the subspace but cannot escape out of it. 
The threshold transfer process [^5) is one example with 
an absorbing subspace. But it does not contain any other 
absorbing state and displays the DP-type transition. 

The structure of absorbing states in the PCPD is 
unique with one point-like absorbing state and one ab- 
sorbing subspace. It is clear that there is no infinite dy- 
namic barrier between these absorbing states. Hence, 
one may argue that the PCPD should belong to the DP 
class. However, this argument turns out to be prema- 
ture. Background diffusing solitary particles generate a 
long-term memory effect on the order parameter (pair 
density). Solitary particles detached from different trains 
of particles diffuse and collide each other to create a new 
particle pair, which leads to history dependence in the 
pair-creation rate. This process is governed by annihilat- 
ing random walks, where the colliding probability of two 
walkers decays algebraically with time. This long-term 
memory effect might be relevant and lead to a new type 
of critical phenomena. 

A numerical investigation using density matrix renor- 
malization group techniques |34j revealed that numerical 
values of some critical exponent ratios are close to the DI 
values rather than the DP values. Subsequent extensive 
simulations |35|, |3(| [37) seem to exclude the possibility 
of both DI and DP classes and suggest that the PCPD 
belongs to a new universality class. However, huge cor- 
rections to scaling conceal the true asymptotic scaling 
behavior and numerical estimates for the critical expo- 
nents are obscure with considerable uncertainty. Simi- 
lar critical behaviors were also observed in related mod- 
els H || |§, which include a modified PCPD |§ with 
the branching process 2X — » 4X instead of 2X —> 3X 
and an effective model with two species of particles cou- 
pled cyclically |38|. 

Despite all those efforts, universal features that could 
characterize the novel universality class were not uncov- 
ered yet. We explore this issue in the present work. Hin- 
richsen noticed that two types of degrees of freedom are 
present in the PCPD; a particle pair that can branch 
and annihilate, and a diffusing solitary particle J38|. All 
activities that can change the number of particles are car- 
ried out by the particle pairs. So the particle pair density 
can be regarded as an order parameter. The two degrees 
of freedom are coupled cyclically; one can be transmuted 
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to the other, and vice versa, via particle diffusion. The 
cyclic coupling results in the long-term memory effect as 
described above, which we believe plays a crucial role in 
this novel critical phenomena. 

These observations lead us to consider a generalized 
PCPD (GPCPD) model with a parameter controlling the 
strength of the memory effect, which will be explained in 
Sec. |l| Without the memory effect, the GPCPD should 
fall into the DP class. So, our parameter connects the 
DP fixed point to the PCPD point, which allows us to 
study the long-term memory effect on the DP-type mod- 
els systematically. 

One may expect that the long-term memory would 
serve as a relevant (at least, marginal) perturbation to 
the DP fixed point, in order to account for the non-DP 
type novel critical behavior at the PCPD point. In that 
case, one can hope that, by varying the control parame- 
ter, those huge corrections to scaling found at the PCPD 
point might be reduced to a tolerable level, so the mod- 
erate numerical efforts may reveal the true asymptotic 
behavior. 

We performed extensive Monte Carlo simulations in 
Sec. [II. Our numerical results show that corrections 



to scaling are huge only near the PCPD point and the 
asymptotic regimes are rather easily reached at other val- 
ues of the control parameter. So, we were able to esti- 
mate the values of the critical exponents with reasonable 
accuracy for a wide range of the control parameter, ex- 
cept not very near the PCPD point. In the absence of 
the memory effect, we find the DP universality class as 
expected. With the memory effect, surprisingly, we ob- 
serve continuous variation of critical exponents with the 
memory strength. Especially, the order parameter expo- 
nent P varies more than 60%, which is far beyond the 
statistical errors less than at most 10%. 



In order to establish this novel universality class with 
continuously varying critical exponents firmly, in Sec. IV, 
we introduce two independent model systems with two 
species of particles: one species plays the role of the par- 
ticle pairs and the other the role of the diffusing soli- 
tary particles. They are coupled cyclically through trans- 
mutations, which leads to the long-term memory effect. 
These models are also studied by Monte Carlo simula- 
tions, and are shown to display the same type of critical 
phenomena as in the GPCPD. 



In Sec. 0, we suggest one possible scenario to account 
for this novel universality class with continuously vary- 
ing exponents. This scenario assumes that the long-term 
memory plays the role of the marginal perturbation to 
the DP fixed point. We numerically check this scenario 
by measuring the life time distribution of solitary parti- 
cles. Our analysis shows that this scenario may be con- 
sidered as a reasonable one, though it needs full field- 
theoretical analysis for a definitive test. Finally, we sum- 



marize our works in Sec. VI 



II. MODEL 

The GPCPD model is defined on a ID lattice of size 
L with the periodic boundary conditions. Each site is 
either occupied by a particle (X) or empty (0). The sys- 
tem evolves in time according to the dynamic rules of the 
PCPD (see Sec. with one additional ingredient: When 
the hopping of a solitary particle creates a new parti- 
cle pair through colliding with another particle, this new 
particle pair annihilates instantaneously with probability 
1 — r or survives with probability r. This does not ap- 
ply to particle pairs formed by branching processes. At 
r = 0, two solitary particles always annihilate upon meet- 
ing and have no chance to turn into a nearest neighbor 
particle pair. The parameter r controls the transmuta- 
tion rate of two solitary particles into a nearest neighbor 
particle pair upon meeting. At r = 1, the model reduces 
to the ordinary PCPD. 

As in the PCP, we take the particle pair density as 
the order parameter. At nonzero r, there is a feedback 
mechanism to increase the pair density via the collision of 
diffusing solitary particles. A solitary particle is created 
from a train of particles by a hopping (detaching) process 
of a boundary particle or a pair annihilation process of a 
triplet of particles. In this process, at least one particle 
pair is sacrificed. The solitary particles diffuse and collide 
with each other to form a particle pair with probability 
r. So, one can say that particle pairs temporarily turn 
into solitary particles and resurrect later by collision of 
those solitary particles. The colliding probability of two 
diffusing particles usually decays algebraically in time. 
Therefore, this feedback mechanism induces long-term 
history dependence, which is called the long-term mem- 
ory effect, of the order parameter. As r increases, the 
long-term memory effect becomes stronger. The PCPD 
model [r = 1) has the maximum memory effect. 

The r = point is special. Collision of two particles 
does not generate a particle pair at all, so there is no 
feedback mechanism for the pair density. Once the sys- 
tem evolves into a state without any particle pair, it stays 
inside a subspace of states without a pair and, in the end, 
will be trapped into one of the absorbing states (the vac- 
uum state and the states with a single particle). So, in 
the viewpoint of the order parameter as the pair density, 
the no-pair subspace serves as the absorbing subspace 
which the system cannot escape from. Evolution inside 
this subspace is governed by the trivial pair annihilation 
dynamics of diffusing particles that have no memory ef- 
fects on the order parameter fluctuations. Therefore, the 
GPCPD at r = should be effectively equivalent to the 
PCP model without diffusion. The no-pair subspace is 
exactly the same as the absorbing subspace of the latter. 
The coupling constants for the branching and annihila- 
tion processes are renormalizcd due to the inclusion of the 
single-particle diffusion, and the critical point is shifted. 
However, the critical fluctuations should be identical and 
we expect the DP-type stationary critical behavior at the 
absorbing transition. We also expect the ordinary DP- 
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FIG. 1: Illustration of dynamic rules. Filled (Empty) circles 
represent occupied (empty) sites. Dots indicate selected sites. 



type scaling in the transient regime, starting from a single 
pair of particles, in contrast to the controversial transient 
behaviors in the PCP model. 

At nonzero r, the system escapes from the no-pair sub- 
space with finite probability. The order parameter can 
be significantly influenced by the memory effect induced 
by this in-and-out dynamics with respect to the no-pair 
subspace. The parameter r connects the DP fixed point 
(r = 0) to the PCPD point (r = 1). Therefore, the 
GPCPD allows us to study systematically the origin of 
the new type of critical behavior found in the PCPD. 

Dynamics of the GPCPD can be implemented in 
Monte Carlo simulations as follows (see Fig. |l|). First, 
select a pair of sites (i, i + 1) at random. When both 
sites are occupied, the two particles (a) annihilate with 
probability (1 — d)p or (b) branch a particle at one of 
the neighboring sites, i — 1 or i + 2, with probability 
(1 — d)(l — p). Branching attempt to an already occu- 
pied site is rejected. When only one site is occupied, (c) 
the particle hops to the other site with probability d. If 
the hopping (not branching) creates a new particle pair, 
the pair (d) survives with probability r or (e) annihilates 
with probability 1 — r. When both sites are empty, no 
change is made. The time increases by one unit after L 
such trials. 

To speed up simulations, we adopted a technique uti- 
lizing a list of active pairs of neighboring sites. A neigh- 
boring site pair is stored in the active pair list if it con- 
tains at least one particle. Then, a site pair in the list is 
selected randomly for dynamics. The time increases by 
1 /Wpajr with iVp a i r the number of active pairs in the list 
at each attempt. 



III. MONTE CARLO SIMULATIONS 

Monte Carlo simulations were performed to investigate 
critical behaviors of the GPCPD. The critical points are 
approached by varying p at each r = 0, 0.25, 0.5, 0.75, 
and 1 with fixed diffusion probability d = 0.1. 



A. Defect simulations 

We performed the so-called defect simulations to locate 
the critical points. Starting with a single pair of parti- 
cles, we measured the survival probability P(t) that the 



system is surviving at time t, the number of particle pairs 
N(t) averaged over all samples, and the mean distance 
of spreading R(t) averaged over the surviving samples. 

Our definition of surviving samples is different from 
the conventional one where all samples not being trapped 
into one of the absorbing states are considered as surviv- 
ing ones. Here, only samples with at least one particle 
pair are regarded as surviving. Samples with only soli- 
tary diffusing particles (no pair but irrespective of the 
number of particles) are regarded as dead ones, even if 
they are not trapped completely. They may be dead for a 
while, but can resurrect later as surviving samples when 
diffusing particles meet and form a particle pair. The 
dead states form the no-pair subspace (see Sec. ||). At 
nonzero r, the system can evolve in and out of the no- 
pair subspace. With this definition, the survival prob- 
ability P(t) actually represents the probability that the 
system contains at least one particle pair, or equivalently 
the probability that the system stays outside the no-pair 
subspace. 

The conventional surviving ensemble includes the most 
states inside the no-pair subspace except the true absorb- 
ing states (the vacuum state and the states with a single 
particle). The survival probability is now dominated by 
the trivial pair annihilation dynamics of diffusing par- 
ticles inside the no-pair subspace [[ll], which does not 
reflect the proper dynamic and stationary scaling behav- 
ior of the order parameter (the pair density). Hence, the 
fluctuations of the order parameter should be described 
in our surviving ensemble defined as the complement of 
the no-pair subspace. 

At criticality, the values of the measured quantities 
decay algebraically as p2| 



P(t) 
N(t) 
R(t) 



t 



-6' 



(1) 



and the double logarithmic plots against time show 
straight lines. Off criticality, these plots show curvatures 
in the long time limit. Precise estimates for the critical 
points and the scaling exponents can be obtained by ex- 
amining the local slopes of the curves. The local slope, 
called as the effective exponent, is defined as 



6'(t) 



Ya[P(t)/P{t/m)) 
lnm 



and similarly for rj{t) and z(t) with a constant m (taken 
to be 10). The power-law scaling behavior implies that 
they converge to the values of <5', 77, and z asymptotically 
as t — > 00 at p = p c . At off-critical points, they behave 
like being at the critical point in the early time regime 
and then deviate to a trivial value characteristic of the 
active or inactive phase. From this crossover behavior 
one can determine the location of the critical point. 

The defect simulations were performed up to 10 5 time 
steps and the observables were measured and averaged 
over ~ 2 x 10 5 samples. Figure || shows plots of the 
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FIG. 2: The effective exponents for the GPCPD at r = 0.5. 



TABLE I: Critical points and critical exponents estimated 
from the defect simulations. As a reference, corresponding 
values of the DP and DI classes, taken from Ref. fea], are 
given. 



r 


Pc 


z 


V 


X 1 









0.04687(2) 


1.58(1) 


0.314(6) 


0.160(5) 


0.160(5) 


0.25 


0.05505(5) 


1.62(3) 


0.29(1) 


0.15(1) 


0.175(5) 


0.50 


0.06636(4) 


1.67(3) 


0.26(1) 


0.14(1) 


0.197(5) 


0.75 


0.08315(5) 


1.75(5) 


0.20(2) 


0.13(2) 


0.235(5) 


1 


0.1112(1) 


1.7(1) 


0.18(5) 


0.09(2) 


0.31 


DP 




1.579(2) 


0.3137(10) 


0.1596(4) 


0.160(2) 


DI 




1.753(3) 


0.000(1) 


0.285(2) 


0.285(4) 



effective exponents against 1/f at r = 0.5. Apparently 
the effective exponents for 5' , 77, and z converge to their 
asymptotic values at p = 0.06636, while they clearly bend 
up or down with a curvature at p = 0.06632 and 0.06640. 
It leads us to estimate that p c = 0.06636(4) and z = 
1.67(3), 77 = 0.26(1), 5' = 0.14(1). 

The errors in the exponent values mainly stem from 
the uncertainty in the p c estimate. Statistical errors are 
much smaller than this systematic error, which can be 
clearly seen in Fig. ||. A correction to the scaling could 
lead to a systematic error in the estimate of p c , and hence, 
of critical exponents. The plots in Fig. || show that the 
correction is quite small for N(t). The effective expo- 
nents S'(t) and z(t) have a little time dependence (seem- 
ingly, linear dependence in 1/f) even at the estimated 
Pc while such time dependence is negligible for 77(f) at 
the estimated p c . So the best estimate for p c is obtained 
from the plot of 77(f) vs. 1/f jl3). We also find similar 
behaviors for other values of r except r = 1. 

We also plot the combination of the effective exponents 
8{t) = l/z(t) - (rj(t) + 5'(t)). At criticality, this expo- 
nent should converge asymptotically to the ratio of two 
stationary scaling exponents, 0/v\\ (known as the initial 
slip exponent), if the hyperscaling relation ]l5| is satis- 
fied. At r = 0.5, we obtain that S = 0.197(5). This value 
will be compared to f}/v\\ estimated independently later 
in the static Monte Carlo simulations. The same analy- 
sis is repeated for other values of r and the results are 
summarized in Table |. 

The two exponents S and 8' do not necessarily coincide 
unless the evolution operator is invariant under the time- 
reversal transformation Q. However, it is well known 
in models with multiple absorbing states that these two 
exponents coincide if one starts with the so-called natural 
initial configurations in defect simulations. This aspect 
will be discussed in details elsewhere [0. 

The PCPD point at r — 1 is an exceptional case. Very 
strong corrections to scaling are observed in all quanti- 
ties as can be seen in Fig. |^. It has been already noted in 
Refs. |}5], ||^]. We tried to locate the critical point from 
the curvature change in the plots of the effective exponent 



vs. 1/t. The plot of 77(f) vs. 1/f shows that p c ~ 0.1113, 
while the plot of 5'{t) vs. 1/f does p c ~ 0.1112. It sug- 
gests that the asymptotic scaling regime has not been 
reached yet until 10 5 time steps at the PCPD point. The 
exponent S shows the worst behavior. Even one can not 
see any difference between super-critical and sub-critical 
behaviors (so we could not estimate an error bar of 6). 
So the values of p c and critical exponents at the PCPD 
limit contain the largest uncertainty. It may be quite in- 
teresting to investigate the origin and the nature of the 
unexpected huge corrections at the PCPD point, but we 
do not have any reasonable explanation at this moment. 

When one faces a strong correction-to-scaling behavior 
in investigating numerically unknown critical phenomena 
of a model system, one may perform an extensive simula- 
tion up to much larger length scales and much longer time 
scales. However, if it is beyond the present day's comput- 
ing capability, one should look for an efficient alternative 
model which presumably exhibits the same critical be- 
havior with smaller correction to scaling. The PCPD 
model applies to such a case. 

The GPCPD model may serve as an efficient alterna- 
tive model for the PCPD model. The r = 1 (PCPD) 
point of the GPCPD model is not a special point; it pos- 
sesses neither an additional symmetry nor conservation 
law, compared to other points for < r < 1. Only the 
strength of the long-term memory effect changes with r. 
One may guess that the GPCPD model with all nonzero 
r may belong to the same universality class as the PCPD 
model and hope that the corrections to scaling are con- 
trollable for small r. It turns out that the GPCPD model 
possesses much smaller corrections to scaling for r < 1 
(at least up to r — 0.75). So the critical points are very 
accurately estimated and the critical exponents are de- 
termined with reasonable errors. 

At r = 0, the values of all exponents are in excellent 
accord with the DP values. It confirms the expectation 
that the GPCPD without the memory effect should be- 
long to the DP class. For other values of r ^ 0, the 
exponent values begin to deviate from the DP values. 
They are also clearly different from the DI values. It 
confirms that the GPCPD with the long-term memory 
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FIG. 3: The effective exponents for the PCPD (r = 1). 



effect (r 0) displays novel type critical phenomena that 
do not belong to the DP nor DI universality class. More 
importantly and very surprisingly, the exponent values 
show a dependence on r (see Table |lj). This opens up a 
new possibility that the novel universality class would be 
characterized by continuously-varying critical exponents, 
i.e., not a fixed point but a fixed line parametrized by the 
memory strength r. This remarkable finding should be 
examined carefully whether it is also present in the scal- 
ing property of the steady states. 



B. Static simulations 

The criticality in the steady states is studied via the 
so-called static Monte Carlo simulations where one starts 
with a macroscopically occupied configuration on finite 
size lattices. Here, we start with the fully occupied con- 
figuration and use the periodic boundary conditions. As 
the order parameter, we measure the pair density p s {t) 
averaged over surviving samples that contain at least one 
particle pair at a given time t, or equivalently averaged 
over the complementary states to the no-pair subspace, 
as defined in the defect simulations previously. The aver- 
aged quantities over the conventional surviving ensemble 
would also describe the system properly in the early time 
regime, up to the time scale when the system typically 
enters into the no-pair subspace. However, since then, 
the trivial pair annihilation dynamics of diffusing parti- 
cles governs the system evolution El| and the average 
value of the order parameter crosses over to the trivial 
stationary value. 

The pair density averaged over our surviving ensemble 
satisfies a finite-size-scaling form as 



(2) 



where L is the system size and e = p c — p is the distance 
from the critical point p c . The exponents and v±_ are 
the critical exponents associated with the order parame- 
ter and the spatial correlation length, respectively. The 
dynamic exponent z is the ratio of the two exponents; 




0.001 



FIG. 4: Left: Order parameter decay at p — p c at each 
r = 0.0 (top),. . . , 1.0 (bottom). Right: The effective expo- 
nent defined with m = 8. Solid lines represent the 
effective exponents at the estimated critical points (see text) 
and dotted (broken) lines at super-critical (sub-critical) points 
at each r = 0.0 (bottom),. . . ,1.0 (top). 



z = v^jvy_ where i/y is the relaxation time exponent (see 
Ref. pi for a review). 

At criticality (p = p c ), we have 



p s (0,t,L) = L-^g(t/L z ) 



(3) 



where g(x) 
for large x. 

Ps ~ t-?/» 

value p s ~ 



~ x~Pl v W for small x and becomes a constant 
Thus, at p — p c , the pair density decays as 
ii for i«L z and saturates to a steady-state 
L-P/vj. f or f>L z . 



The power-law scaling behavior in the transient 
regime (i <C L z ) can be studied with an effective ex- 
ponent 

io, i_ HPs(t)/p s (t/m)} 

\P/ V \\\ = ; 

mm 

with a constant m. It converges to ft/vn for large t (<C 
L z ) at the critical point p = p Cl but deviates from it 
at p ^ p c . Using this property we could determine the 
critical point p c independently and hence the exponent 

Pin- 

The order parameter is measured in a lattice of size 
L = 10 4 up to t = 10 5 time steps and averaged over 
2000 ~ 5000 samples. The finite-size saturation effect 
is invisible up to t = 10 5 for this system size. Fig- 
ure ^ shows the log-log plots of the order parameter 
at the estimated critical points as well as the plots of 
the effective exponents at and near the criticality. From 
these plots, we estimate p c = 0.04687(2), 0.05505(5), 
0.06637(2), 0.08317(4), and 0.1113(1) at each value of 
r = 0, 0.25, 0.5, 0.75, and 1, respectively. Except for 
the r = 1 case, the values of the effective exponents at 
the estimated p c do not show any significant systematic 
error and seem to reach their asymptotic values quite 
early, at around t = 10 3 . Statistical fluctuations are also 
very small. The off-critical data in Fig. |4| at the values 
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TABLE II: Critical exponents of the GPCPD and correspond- 
ing values of the DP and DI classes, taken from Ref. p2|. 



r 


HI "|| 


(3/v ± 


z 









0.159(1) 


0.252(3) 


1.58(1) 


1.10(1) 


0.277(5) 


0.25 


0.173(5) 


0.282(5) 


1.64(5) 


1.10(3) 


0.310(14) 


0.50 


0.197(3) 


0.330(6) 


1.69(5) 


1.10(3) 


0.363(17) 


0.75 


0.230(5) 


0.40(1) 


1.72(5) 


1.17(5) 


0.468(30) 


1 


0.27(4) 


0.50(4) 


1.8(2) 


1.30(10) 


0.65(12) 


DP 


0.1596(4) 


0.2522(6) 


1.5798(18) 


1.0972(6) 


0.2767(4) 


DI 


0.285(2) 


0.500(5) 


1.750(5) 


1.84(6) 


0.92(3) 
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of p displaced by the amount of error bars in the above 
p c estimates, clearly show either upward or downward 
curvatures, which demonstrates the accuracy of our p c 
estimate. Moreover, these p c values are in very good 
agreement with those obtained from the defect simula- 
tions (see Table |). 

The values of P/v\\ are given by the limiting values of 
the effective exponents for large t at the estimated p c , 
which are presented in Table [II]. For consistency, we use 
the defect simulation results for the p c values and their 
errors. The estimated values of &fv\\ are again in excel- 
lent accord with the values of S — l/z—(r) + 6') measured 
in the defect simulations for 1, which implies that the 
hyperscaling relation jl5| holds in this model. 

Note that the correction to scaling is significant at the 
PCPD point (r = 1), as found in the defect simulations. 
At r = 1 , the effective exponent plot shows a strong time 
dependence even at the estimated critical point and there 
is a noticeable curvature in Fig. [|. It leads to a rather 
large error in the p c estimate and, hence, the exponent 
estimates. Our estimate of also shows a noticeable 
discrepancy with the estimate of 5 (see Tables | and ||) . 
However, this discrepancy presumably comes from the 
inaccurate estimation of p c at r — 1 and insufficient time 
steps getting into the asymptotic regime. So it cannot 
be regarded as an evidence for the breakdown of the hy- 
perscaling relation at the PCPD point. 

The steady-state pair density for t 3> L z satisfies 



p s {e, t = oo, L) = L 



-P/"± 



h(eL 



(4) 



where the scaling function behaves as h(x) ~ x@ for large 
x and becomes a constant near x ~ 0. At criticality 
(e = 0), it decays algebraically with size as p s ~ L~P' VJ - . 
Utilizing this algebraic scaling property, we could also 
estimate the location of the critical points p c . They again 
agree perfectly well with the previous other results. 

We run 10 3 ~ 10 4 samples up to 4 x 10 6 time steps 
for the system size L = 2 6 , . . . , 2 11 at the estimated crit- 
ical points given in Table ffl. Plots of p s vs. L at each 
critical point are presented in Fig. |^. A least-square 
straight-line fitting in the log-log plot might yield the 
value of f3/is±. However, it would result in an inaccu- 
rate estimate in the presence of a strong correction to 



FIG. 5: The power-law scaling behavior of the order param- 
eter at p = p c at each r. The solid line is a guide for the eyes 
whose slope is @/v±_ in Table [o]. 




t I L 

FIG. 6: Scaling plots according to Eq. (B) at each r = 0.0 
(top),. . . , 1.0 (bottom). Each data set is shifted vertically by 
a constant factor to avoid an overlap with others at different 
values of r. 



scaling. Instead, we made use of an effective exponent 
[/3/V±] = ~ m [p s (2i)//5 s (i)]/ m 2 to extract the accu- 
rate exponent value from the extrapolation. It is plotted 
in the inset of Fig. || and the results are presented in Ta- 
ble [n| As before, there exists a strong correction to the 
power-law scaling at the PCPD point at r = 1, whereas 
such an effect is very small for the other values of r. 

With p c and 0/v± determined, we could estimate the 
dynamic exponent z using the scaling form in Eq. (^|). It 
was estimated as the optimal value that yields the best 
collapse of p s (s = 0, t, L) for L = 2 7 , . . . , 2 11 (see Fig. |). 
The resulting values of z are presented in Table |l[ They 
are consistent with the ratio of the two exponent ratios, 
P/v± and (3/v\\ , and agree very well with the values mea- 
sured in the defect simulations. 

We also determined the value of v± by collapsing off- 
critical steady-state data of p s (e, t = oo, L) in the scaling 
plot, using Eq. ([I|). In Fig. ^|, the data for each r collapse 
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FIG. 7: Scaling plots according to Eq. (^) at each r — 0.0 
(bottom),. . . , 1.0 (top). The solid line is a guide for the eyes 
whose slope is (3 in Table [o]. Each data set is shifted verti- 
cally by a constant factor to avoid an overlap with others at 
different values of r. 



well on a single curve with a specific value of v±_ , respec- 
tively. The order parameter exponent (5 is then obtained 
by a product of j3/v±_ and v±_ and presented in Table 0. 
They are equal to the slopes of the scaling function h(x) 
for large x in the log-log scale as shown in Fig. [?]. 

At the PCPD point (r = 1), our results are consistent 
with those obtained by Carlon et al. Q (p c = 0.111(2), 
z = 1.87(3), 0/v± = 0.50(3)) and Hinrichsen (35) (p c = 
0.1112(1), /3/z/|| = 0.25(2), z = 1.83(5), (3 < 0.67). Al- 
though some of critical exponent ratios are close to the 
DI values, the critical exponents are incompatible with 
those of the DI class. They are also inconsistent with the 
DP universality class. The strong correction-to-scaling 
behaviors were observed at the PCPD point with r = 1 
in both the defect and the static simulations. So our esti- 
mates for the critical exponents contain considerable un- 
certainty at the PCPD point. However, such corrections 
are not prominent for other values of r, which enables us 
to estimate the values of the critical exponents very ac- 
curately for r < 1. As expected, the values of the critical 
exponents at r — agree very well with the DP values. 

For < r < 1, the correlation length exponent i/± 
seems to remain the same as the DP value for a wide 
range of r, at least up to 0.75. The value of the relaxation 
time exponent i/u = zis± seems to increase slightly with 
increasing r (10% up to r — 0.75), but, within present 
numerical accuracy, it may be difficult to conclude that 
its variation is real and not due to corrections to scaling. 
However, it is clearly visible that the order parameter ex- 
ponent (3 (also 0/v± and varies considerably with 
r (more than 60% !). It clearly signals a novel univer- 
sality class with scaling exponents varying continuously 
with the parameter r. 

Summing up all results, our Monte Carlo simulations 
show that the GPCPD displays the novel critical phe- 
nomena distinct from the DP and DI classes, with con- 
tinuously varying exponents depending on the strength 



of the long-term memory effect r. These results are quite 
surprising and remarkable. In the renormalization group 
language, it implies that there is a fixed line parametrized 
by r, instead of fixed points at two end points; the DP 
(r = 0) and the PCPD point (r = 1). 

Continuously varying critical exponents in non- 
equilibrium systems are very rare. Systems with in- 
finitely many absorbing states like the PCP may dis- 
play continuously varying exponents depending on ini- 
tial conditions, but only the exponents describing non- 
stationary properties (<$' and rj) are varying Jh3[. This 
variance is even disputed very recently jlT], |18| . Levy- 
type long-range flights are known to be relevant to ab- 
sorbing critical phenomena and the stationary critical 
exponents vary continuously with the exponent describ- 
ing the long-range tail distribution of the flights [^7| . 
However, the GPCPD involves only short-range pro- 
cesses. Recently, multi-species particle reaction-diffusion 
systems with inter-species hardcore interactions are con- 
jectured to display continuously varying stationary-state 
exponents, but the transition occurs at the trivial anni- 
hilation point ^8|. Moreover, the variance is very small 
and it would be extremely difficult to confirm it numer- 
ically. Our studies suggest strongly that the GPCPD 
belongs to a completely new type of universality class 
with continuously varying exponents which has not been 
explored before. 



IV. UNIVERSALITY 

In order to establish the novel universality class firmly, 
we study other model systems that share the common 
feature of the memory effect with the GPCPD. As de- 
scribed in Ref. p8|, the PCPD can be regarded in a 
coarse-grained level as a cyclically coupled system of two 
particle species: one species (A) performing DP-like dy- 
namics and the other species (B) the annihilating ran- 
dom walks. In the GPCPD model, a particle pair cor- 
responds to an A particle and a solitary particle to a B 
particle. Two species are coupled through transmuta- 
tions, which leads to the long-term memory effect. We 
set up two different ID models of a cyclically coupled 
system and investigate numerically their scaling behav- 
ior to check our novel universality class. The two typical 
DP dynamics are employed for the A-particle dynamics, 
i.e. the branching-annihilating random walk model with 
one offspring (BAW1) |l2) and the contact process. 

The first model (referred to as the ABB model) adopts 
the BAW1 evolution rule for the A dynamics. Each lat- 
tice site is either occupied by an A or B particle, or 
empty. First, select a site at random. When it is oc- 
cupied with an A, the A particle (a) branches an A with 
probability (1 — p) on a neighboring site, or (b) hops to 
a neighboring site with probability p(l — fi), or (c) trans- 
mutes to a B particle with probability pfi. When the 
branched or hopping A particle lands on another A par- 
ticle, both A particles annihilate immediately as in the 
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BAW models. In case that it would land on a B par- 
ticle, the trial is rejected and the B particle transmutes 
to the A particle. When the selected site is occupied 
with a B particle, the B particle hops to a neighboring 
site with probability d or does nothing with probability 
1 — d. When it would land on another particle (cither 
A or B), the trial is rejected and all involved B parti- 
cles transmute to A particles. When the selected site is 
vacant, nothing happens. 

The evolution rule of the second model (the ABC 
model) is almost equivalent to the ABB model, except 
that the contact process is adopted for the A particle 
dynamics. An A particle (a) branches an A with proba- 
bility (1 — p) on a neighboring site, or (b) vanishes spon- 
taneously with probability p{\ — p), or (c) transmutes to 
a B particle with probability pp. When the branched A 
particle would land on cither A or B particle, the trial 
is rejected and the involved B particle transmutes to the 
A particle. A B particle hops to a neighboring site with 
probability d or does nothing with probability 1 — d, as in 
the ABB model. When it would land on another particle 
(either A or B), the trial is rejected and all involved B 
particles transmute to A particles. 

In both models, the two species of particles are cou- 
pled cyclically through the transmutations. The trans- 
mutation rate from A to B is controlled explicitly by the 
parameter (i. As p increases, the transmutation events 
occur more likely . At p = 0, the A^> B channel is com- 
pletely blocked. Therefore there is no feedback mecha- 
nism to change the A particle density through interme- 
diate B particles. The transmutation rate from B to 
A is implicit and determined by the evolution rule and 
the other control parameters. The B — > A channel re- 
lies heavily on the the diffusive property of B particles. 
With nonzero d, the B particles diffuse until they meet 
another particle and transmute to the A particles. This 
process invokes the same mechanism that gives rise to the 
long-term memory effect in the GPCPD model for /i ^ 0. 
The B particles created by the A particles through trans- 
mutations perform the random walks before transmuting 
back to the A particles at later times. At d = 0, the 
B — > A channel is still open, but only with the short- 
term memory processes, which should be irrelevant. 

We expect to observe the same type of critical behav- 
iors, as in the GPCPD model, characterized by the con- 
tinuously varying critical exponents. We performed the 
defect and static simulations to locate the critical points 
and estimate the critical exponents of the ABB and the 
ABC models at several values of (i with fixed hopping 
probability d = 1. In Fig. ||, we compare the space- 
time structures of spreading patterns of activities in the 
GPCPD (r = 0.5), ABB (p = 0.5), and ABC (p = 0.2) 
models at critical points. In all cases, the intermedi- 
ate long-range diffusive processes of the solitary particles 
(or the B particles) are commonly observed between the 
particle pairs (or A particles) at different space-time po- 
sition. The time scale of these processes represents the 
effective life time of the intermediate B particles, that 
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FIG. 8: Critical spreading of activities originated from a 
seed (a particle pair for the GPCPD and an A particle for 
the ABB and ABC models). A black pixel represents a par- 
ticle pair or an A particle, and a gray one a solitary particle 
or a B particle. The horizontal (vertical (down)) direction 
corresponds to the space (time) direction. 
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FIG. 9: Parametric plots of /3/v± (filled symbol) and (t/v\\ 
(open symbol) with respect to (3 of the GPCPD, the ABB, 
and the ABC models. The data taken from Refs. |36[ |3$], 
p9| , ffo| are also plotted and labelled by "Others" . 



seems to be comparable to the simulation time. This as- 
pect will be discussed quantitatively in the next section. 

We take the A particle density as the order parameter 
and use the same definition of the surviving samples as in 
the GPCPD model, i.e. only the samples with at least one 
A particle are regarded as surviving ones. Here, we only 
state the results without showing the data. Our numeri- 
cal simulations confirm that both ABB and ABC models 
show continuously varying exponents with the parameter 
p. As expected, both models converge to the DP class at 
(i = 0. 

In Fig. |[ the values of the critical exponents at vari- 
ous transmutation rates are plotted together with those 
for the GPCPD model and those available from previous 
other works |35|, [36], [38| [h| . Remarkably, they are ly- 
ing along a single smooth line. It indicates that all these 
models fall into the same novel universality class with 
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FIG. 10: Log-binned life time distribution of solitary particles 
at the critical points of the GPCPD with r = 0.25, . . . , 1.00. 
Inset: Effective exponent for F(r) ~ t~° . The apparent cur- 
vature for large r > 10 3 is simply due to the finite simulation 
time of i m ax = 10 4 . 



continuously varying exponents that can be parametrized 
by a single parameter (memory strength). The lines for 
/3/v± and P/vn seem to be almost linear for small (3 < 0.5. 
It implies that the values of the correlation exponents, vj_ 
and v\\ , do not vary too much from the DP values, in con- 
trast to a wide variation of the order parameter exponent 
(3. 



V. LIFE TIME DISTRIBUTION OF 
INTERMEDIATE PARTICLES 

The continuously varying exponents of single parame- 
ter family remind us of the equilibrium Gaussian model 
in two dimensions where a marginal operator accounts for 
the continuous variation of the critical exponents. Here, 
we test a similar scenario that the long-term memory ef- 
fect mediated by the solitary particles (B) may play the 
role of the marginal operator. 

In a heuristic way, one can write down a mean field 
type rate equation for the order parameter, the pair (A) 
density pit), as 

^ = ap(t) - bp{tf + A J* dt'F(t - t')p(t>) , (5) 

where the first and the second terms appear generically 
representing the spontaneous annihilation and branching 
of the pairs (A). The A term represents the transmuta- 
tion of the solitary particles (B), which were detached 
(created) from A at time t', into A at later time t > t'. 
Therefore, the kernel F(t) represents the distribution 
function of the life time r of the B particles. 

We expect that F(t) should have a power-law type 
distribution, because the transmutation event from B to 
A is governed by an annihilating random walk process 
of two B particles. We measure numerically F(t) at the 



critical point from the histogram of life time r of solitary 
(B) particles. We start with a single A on an infinite 
lattice, and generate 10 5 samples up to i max = 10 4 to 
obtain at least 10 9 B's. Figure |l^ shows F(t) averaged 
over surviving samples for the GPCPD model. It follows 
a power law 

F(r) ~ r- e . (6) 

The values of 9 are determined from the effective expo- 
nent plots as shown in the inset. After the transient 
regime at small t, the effective exponent converges to 
= 2.25(5) for all values of r. 

The apparent deviation from the power-law distribu- 
tion for large r > 10 3 is simply due to the finite simula- 
tion time of £ m ax = 10 4 . The number of B particles with 
a long lifetime comparable to i max must have been un- 
derestimated in our simulations, which leads to a slight 
downward deviation of F(t) for large r. Therefore, the 
apparent blowup of the effective exponent 9 near 1/r — 
in the inset of Fig. [l0| should be ignored. The power law 
distribution with 9 = 2.25(5) is observed universally for 
all three models considered in this paper. Note that the 
power-law distribution sets in later as one approaches 
the PCPD point at r = 1. Presumably this is related to 
the strong correction-to-scaling behavior observed in the 
PCPD. 

A straightforward power counting analysis of the mean 
field equation (^) shows that the renormalization group 
(RG) scaling exponent of the A term is given by y\ — 
z(2 — 9) with z the dynamic exponent. With 9 > 2, 
the long-term memory effect appears to be irrelevant 
and our scenario for the marginal operator may be dis- 
carded. However, this naive mean-field-typc analysis may 
go wrong in the interacting theory with fluctuations con- 
sidered. The Levy-flight DP system is such an exam- 
ple p7| , where it is analytically shown by the e-expansion 
type RG analysis that the relevancy of long-range flights 
sets in a little bit earlier than expected from the mean 
field theory. We suspect that a similar situation also oc- 
curs in the GPCPD and the marginal value of 9 becomes 
slightly bigger than 2. In this point of view, our sce- 
nario may be still alive and even considered reasonable 
because 9 = 2.25(5) is not far away from 2. A definitive 
test on our scenario needs accurate determination of the 
marginal value in a full field-theoretical context, which is 
beyond the scope of the current paper. 

VI. SUMMARY AND DISCUSSION 

In summary, we introduced a generalized version of 
the PCPD (GPCPD) with a parameter controlling the 
long-term memory effect. The GPCPD connects the DP 
fixed point to the PCPD point continuously. We inves- 
tigated numerically the nature of the absorbing phase 
transitions for the GPCPD in one dimension. Our nu- 
merical results strongly suggest that the GPCPD belongs 
to the novel universality class, which is characterized by 
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the long-term memory effect and the continuously vary- 
ing critical exponents. This model can be viewed as the 
cyclically coupled systems based on the DP dynamics 
with the long-term memory effect. We showed numeri- 
cally that the two other variants of these systems fall into 
the same universality class. 

We presented one possible scenario to account for this 
novel universality class that the intermediate particles 
generate the long-term memory effect which may play 
the role of the marginal operator to the DP fixed point. 
We measured the lifetime distribution of the intermedi- 
ate particles, which turned out to follow the power law 
with the exponent value slightly larger than the marginal 
value in a naive noninteracting theory. This small differ- 
ence may disappear in a full interacting theory and the 
long-term memory effect may become really marginal. 
However, this is only speculative and definitely needs a 
full field-theoretical treatment. 

Very recently, there have appeared several works on the 
PCPD and related models. Park and Kim Jl5| studied 
three different variants of the PCPD along a special line 
in the parameter space of the diffusion and the reaction 
rates. They found numerically that the exponent values 
for all three models are consistent within statistical er- 
rors, e.g., P/v\\ ~ 0.24(1). However, this does not imply 
that the PCPD scaling with general parameter values can 
be also described by the same fixed point. 

Dickman and de Menezes [fl6| studied the PCPD using 
the same ensemble for the surviving samples as we used 
here. They found that this ensemble (called as the reac- 
tive sector) is quite useful to study the finite size scaling 
of the order parameter. They observed numerically that 
the critical exponents vary with the diffusion rate d. Sim- 
ilar to our results, the correlation exponent v± appears 
to be independent of d and the same as the DP value. 
The relaxation time exponent v\\ and the order parame- 
ter exponent (3 vary with d up to only about 20% (much 
smaller than the variation in the GPCPD ~ 60%), which 
leads to an indecisive conclusion. 

More recently, Hinrichsen |47| introduced a cellular au- 
tomaton, which presumably belongs to the same uni- 
versality class as the PCPD. Using a parallel update, 
he simulated the model with size up to L — 2 21 until 
t = 2.5 x 10 6 . This work confirms again the difficulty 
in studying the critical behavior of the PCPD. The sys- 
tem does not reach the asymptotic scaling regime at that 
time scale. From the temporal trend of the effective ex- 
ponents, he suggested an extremely slow crossover to DP, 
but this suggestion is also far from being conclusive. 

Kockelkoren and Chate [Q studied general reaction- 
diffusion processes without the fermionic constraint. 
From Monte Carlo study up to L = 2 22 and t ~ 10 7 , they 
obtained — 0.200 in a model characterized with 

XX -> XXX and XX -» 0. They reported that the 
same exponents are observed in other similar models, and 
hence claimed that there is a single universality class for 
the PCPD. However, their study does not cover the gen- 
eral cases. Odor EtI studied the PCPD with Monte Carlo 



simulations up to L = 10 5 and t = 10 s , and obtained 
/3/v\\ — 0.21 for high diffusion rates and — 0.25 fol- 
low diffusion rates. On the other hand, assuming a log- 
arithmic correction, Odor obtained (3/v\\ — 0.21 also for 
the low diffusion rates. 

These numerical results seem to favor a single non-DP 
universality class for the PCPD. However, as hinted in 
the work of Dickman and de Menezes Gq], there may be, 
if any, a rather small variation (< 20%)of the scaling ex- 
ponents by changing diffusion rates. So it is not surpris- 
ing to see an apparent single universality under the effect 
of huge corrections to scaling. In contrast, the exponent 
variations in the GPCPD are much bigger (~ 60%), so it 
is easier to confirm their variations. 

We believe that the long-term memory effect is respon- 
sible for unusually long relaxation and non-DP scaling, 
and possibly for continuously varying exponents. The 
long-term memory is also present in the ordinary PCPD 
model, but it is controlled implicitly via diffusion rates 
d of isolated or intermediate particles. So, it may not 
be clear to predict how the long-term memory emerges 
with d. In the GPCPD model, we directly control the 
strength of the long-term memory effect by varying the 
parameter r. 

The GPCPD model smoothly connects the DP fixed 
point (r = 0) to the PCPD model (r — 1), in a sense 
that the exponent values change with r continuously and 
monotonically, starting from the DP values. In contrast, 
the PCPD exponent values seem to jump from the DP 
values and decrease slightly with the diffusion rate d. 
For example, as soon as we turn on the diffusion process, 
/3/v\\ increases abruptly from the DP value of 0.1595 to 
~ 0.25 at d = 0.1 and slightly decreases to 0.20 ~ 0.23 
at high d = 0.7 ~ 0.8 |f || §|. Even if one accepts 
the claim that this variation is not real and /3/^|| con- 
verges to ~ 0.20 asymptotically at any nonzero d (single 
universality class), one can not avoid the fact that cor- 
rections to scaling are much bigger at low d than at high 
d. This suggests that the d — > limit in the PCPD 
model should encounter an unusual crossover behavior. 
In the ordinary crossover, one expects an interference of 
the DP fixed point on the PCPD universality class at 
d ~ 0, which generates an apparent exponent value be- 
tween the DP and the PCPD value. All numerical results 
simply disagree with this. The apparent value (0.25) for 
(3/v\\ at low d is higher than both the DP (0.16) and the 
estimated PCPD value (0.20). Therefore, the d = limit 
in the PCPD model is unusually singular, which makes it 
impossible to study the PCPD scaling behavior system- 
atically starting from the DP fixed point. In contrast, 
our GPCPD model is generically well suited to a system- 
atic investigation of the PCPD model by controlling the 
long-term memory effect directly. 

Besides that, there is one important technical point we 
adopted in this paper for analyzing the numerical data. 
As explained in Sec. Ill, we chose the surviving ensem- 



ble as the collection of samples with at least one par- 
ticle pair and the order parameter as the particle pair 
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density. With the conventional choice for the surviving 
ensemble (samples not trapped in one of the absorbing 
states), we found that the order parameter in finite size 
systems bears two time scales (nontrivial relaxation time 
and trivial pair annihilation time) and does not show 
simple scaling collapse with one- variable scaling function 
like in Eq. (^|) p4| . Our choice of the surviving sample 
should involve only one time scale and leads to an excel- 
lent scaling as shown in FigjFj. With this ensemble, it is 
natural to choose the particle pair density (rather than 
the particle density) as the order parameter. 

Summing up, we believe that the GPCPD model serves 
as an efficient generalized model to study the PCPD 
model systematically. The GPCPD seems to have much 
less corrections to scaling, compared to the ordinary 
PCPD model, which enables us to present accurate nu- 
merical data for the exponent values. In Figs. || and ||, 
one can clearly see that the asymptotic values for the ex- 
ponents set in quite early, at around t = 10 3 , except the 



PCPD case (r = 1). In order to check any possible long 
time crossover, we performed an extra static simulation 
for the GPCPD at r = 0.25 on a lattice of size L = 10 5 up 
to t = 10 7 . We found the more accurate critical point at 
p c = 0.055045(5) and the exponent j3/v\\ ~ 0.172(2) g| . 
The result is in excellent accord with our early-time resu t 
on a smaller lattice (see Table | and [n|). This confirms 
again that the GPCPD at r ^ 1 does not suffer from 
huge corrections to scaling as observed in Refs. |47], |49| . 



Acknowledgments 

We thank P. Grassberger for his critical reading of this 
manuscript and valuable comments. We also thank Su- 
Chan Park for useful discussions. This work was sup- 
ported by grant No. 2000-2-11200-002-3 from the Basic 
Research Program of KOSEF. 



[1] For a review, see J. Marro and R. Dickman, Nonequi- 

librium Phase Transitions in Lattice Models (Cambridge 

University Press, Cambridge, 1999). 
[2] H. Hinrichsen, Adv. Phys. 49, 815 (2000). 
[3] P. Grassberger, F. Krause, and T. von der Twer, J. Phys. 

A 17, L105 (1984); P. Grassberger, J. Phys. A 22, L1103 

(1989). 

[4] H. Park and H. Park, Physica A 221, 97 (1995). 

[5] H. Hinrichsen, Phys. Rev. E 55, 219 (1997). 

[6] W. Hwang, S. Kwon, H. Park, and H. Park, Phys. Rev. 

E 57, 6438 (1998). 
[7] H.K. Janssen, Z. Phys. B 42, 151 (1981). 
[8] P. Grassberger, Z. Phys. B 47, 365 (1982). 
[9] T.E. Harris, Ann. Prob. 2, 969 (1974). 
[10] E. Domany and W. Kinzel, Phys. Rev. Lett. 53, 311 
(1984). 

[11] R.M. Ziff, E. Gulari, and Y. Barshad, Phys. Rev. Lett. 

56, 2553 (1986). 
[12] H. Takayasu and A.Yu. Tretyakov, Phys. Rev. Lett. 68, 

3060 (1992). 

[13] I. Jensen, Phys. Rev. Lett. 70, 1465 (1993); I. Jensen and 
R. Dickman, Phys. Rev. E 48, 1710 (1993). 

[14] G. Odor, J.F. Mendes, M.A. Santos, and M.C. Marques, 
Phys. Rev. E 58, 7020 (1998). 

[15] J.F.F. Mendes, R. Dickman, M. Henkel, and M.C. Mar- 
ques, J. Phys. A 27, 3019 (1994). 

[16] R. Dickman and D. ben-Avraham, Phys. Rev. E 64, 
020102 (2001). 

[17] P. Grassberger, H. Chate, and G. Rousseau, Phys. Rev. 

E 55, 2488 (1997). 
[18] A. Jimenez-Dalmaroni and H. Hinrichsen (unpublished). 
[19] F. van Wijland, Phys. Rev. Lett. 89, 190602 (2002). 
[20] N. Menyhard, J. Phys. A 27, 6139 (1994); N. Menyhard 

and G. Odor, ibid. 28, 4505 (1995). 
[21] M.H. Kim and H. Park, Phys. Rev. Lett. 73, 2579 (1994); 

H. Park, M.H. Kim, H. Park, Phys. Rev. E 52, 5664 

(1995). 

[22] I. Jensen, Phys. Rev. E 50, 3623 (1994). 



[23] N. Inui and A. Yu. Tretyakov, Phys. Rev. Lett. 80, 5148 
(1998). 

[24] W. Hwang and H. Park, Phys. Rev. E 59, 4683 (1999). 
[25] M.C. Marques and J.F.F. Mendes, Eur. Phys. J. B 12, 
123 (1999). 

[26] H. S. Park and H. Park, J. Korean Phys. Soc. 38, 494 

(2001) . 

[27] H.K. Janssen, K. Oerding, F. van Wijland, and H.J. Hil- 
horst, Eur. Phys. B 7, 137 (1999); H. Hinrichsen and M. 
Howard, ibid. 7, 635 (1999). 

[28] S. Kwon, J. Lee, and H. Park, Phys. Rev. Lett. 85, 1682 
(2000). 

[29] M. Rossi, R. Pastor-Satorras, A. Vespignani Phys. Rev. 

Lett. 85, 1803 (2000). 
[30] A. Lipowski and M. Droz, Phys. Rev. E 65, 056114 

(2002) . 

[31] A.J. Noest, Phys. Rev. Lett. 57, 90 (1986); M. Bramson, 
R. Durrett, and R. H. Schonmann, Ann. Prob. 19, 960 
(1991); A.G. Moreira and R. Dickman, Phys. Rev. E 54, 
R3090 (1996) ; R. Cafiero, A. Gabrielli, and M. Munoz, 
Phys. Rev. E 57, 50 60 (1998); J. Hooyb erghs, F. Igloi, 
and C. Vanderzande, |sond-mat/020361C| . 

[32] A binary spreading process like the PCPD was first in- 
troduced and studied by Grassberger in ||. 

[33] M. J. Howard and U. C. Tauber, J. Phys. A 30, 7721 
(1997). 

[34] E. Carlon, M. Henkel, U. Schollwock, Phys. Rev. E 63, 

036101 (2001). 
[35] H. Hinrichsen, Phys. Rev. E 63, 036102 (2001). 
[36] G. Odor, Phys. Rev. E 62, R3027 (2000). 
[37] H. Chate and P. Grassberger (unpublished). 
[38] H. Hinrichsen, Physica A 291, 275 (2001). 
[39] K. Park, H. Hinrichsen, and I.-M. Kim, Phys. Rev. E 63, 

065103 (2001). 
[40] G. Odor, Phys. Rev. E 63, 067104 (2001). 
[41] This trivial dynamics may be "dressed" by the intefer- 

ence with the local clouds of particle pairs. This issue will 

be discussed elsewhere. 



13 



[42] P. Grassberger and A. De La Torre, Ann. Phys. 122, 373 [45] 
(1979). [46] 

[43] It does not imply that the estimate using the plot of S'(t) 

or z(t) vs. 1/t would yield a different value of p c . In such [47] 
cases the upward and downward curvature at offcritical [48] 
points guides the estimation, which leads to the same [49] 
result. 

[44] J. D. Noh and H. Park (unpublished). 



K. Park and I.-M. Kim, Phys. Rev. E 66, 027106 (2002). 
R. Dickman, M.A.F. de Menezes, Phys. Rev. E 66, 
045101(R) (2002), 



H. Hinrichsen, |cond-mat/0208345| (2002) 



J. Kockelkoren and H. Chate 
G. Odor 



cond-mat/0209287 (2002) 



cond-mat/0208497 (2002) 



